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Abstract 

The thermal contact conductance of a 0.4 um 
surface finish OFHC copper sample pair has been 
investigated from 1.6 to 3.8 K for a range of 
applied contact forces up to 670 N. 

Experimental data have been fitted to the 
relation 

Q = J*aT n dT 

by assuming that the thermal contact conductance 
is a simple power function of the sample 
temperature. 

It has been found that the conductance is pro- 
portional to T 2 and that conductance increases 
with an increase in applied contact force. These 
results confirm earlier work. 


Introduction 

The optimum design of cryogenic instruments 
requires accurate thermal models. This is espe- 
cially important for instruments where performance 
is sensitive to temperature. Infrared instruments 
such as the Infrared Astronomical Satellite (IRAS) 
and the Shuttle Infrared Telescope Facility 
(SIRTF) fall into this category. The present 
models are limited by a lack of knowledge of the 
low temperature thermal conductance of the bolted 
joints that are typically used in the instrument- 
to-system interface. Previous studies of pressed 
contacts though limited in scope have shown that 
the thermal conductance does not obey the 
Wiedemann-Franz law (that states that the ratio of 
thermal to electrical conductivities is propor- 
tional to the temperature). In this paper, an 
effort to characterize the thermal conductance of 
pressed contacts at liquid helium-4 temperatures 
is described. Specifically, the dependence of 
thermal contact conductance on applied force and 
temperature is discussed. 


Theory and Previous Work 

The phenomenon of thermal contact resistance 
is attributable to several factors; most notably, 
it is the consequence of contact being made only 
at discrete locations, rather than over the entire 
surface area. Ideally, the contact area is repre- 
sented by the interface area of the surfaces them- 
selves; however, a close examination reveals that 
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even in the smoothest surfaces, irregularities 
exist which may restrict the contact area to as 
few as three discrete spots, irrespective of the 
sample dimensions. This theory is supported by 
experimental findings that contact resistance is 
dependent on the force of contact rather than con- 
tact area. 1 * 2 As the contact force increases, the 
material deforms. Thus, the area of initial con- 
tact increases and new spots develop. The heat 
flow is constricted in the vicinity of the contact 
locations because of the narrowness of the effec- 
tive areas of contact. This constriction is in 
large part responsible for the contact resis- 
tance. 3 * 4 Estimates of the constriction resistance 
have been made for various assumed contact geom- 
etries by modeling the contacts as individual 
elements. 5 By arranging the elements in groups of 
varying heights, the case of surface waviness can 
be accounted for as well. 3 ’ 5 

A relation involving contact pressure and 
material hardness exists for determining the ratio 
of the surface area to the actual area of contact. 
However, the equivalent radius of the contact spot 
must be known and, at the present, there is no 
theoretical method available for general determina- 
tion since each sample must be considered on an 
individual basis. 4 The height of the contact gap 
is also significant and a method does exist for 
its estimation. 3 

Additionally, the presence of surface films 
or oxides contribute to the problem. This is 
especially significant in the case of materials 
such as aluminum, which form oxides immediately in 
the presence of air. In the case of oxides, the 
oxide layer must be penetrated to obtain a consis- 
tent measure of the thermal resistance. 

Although there has been significant inter- 
est 3,5 ~ 10 in the problem of contact resistance, 
estimation of the resistance from proposed theo- 
retical models is still not a simple task. Princi- 
ples of variational calculus have been applied to 
the problem to determine upper and lower bounds; 
however, because of the nature of the problem, 
there are limitations in the models and no satis- 
factory agreement exists between models due to the 
incompatibility of their boundary conditions. 

This poses a problem in predicting the behavior of 
pressed contacts and, therefore, most usable data 
in the field are empirical. Additionally, an 
attempt to predict thermal conductance from the 
electrical conductance given by the Wiedemann- 
Franz relation yields values much less than those 
obtained empirically. It has been found that this 
discrepancy can be as large as a factor of 
5 x 10 6 for contacts made at 4.2 K (Ref. 1). 

Previous work consists of limited data for 
pressed contacts in the 4 K region. 1 ,2 * 11-14 Of 
these data, Cu-Cu results seem to dominate and, of 
the work surveyed, there seems to be good general 
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agreement for the Cu-Cu conductance under applied 
forces of 450 N. Other work is limited to particu- 
lar sample pairs and configurations, often corre- 
sponding to specific applications. 9 * 15 A need 
exists for more general thermal conductance data 
covering a variety of samples over a range of tem- 
peratures, contact forces, and surface conditions. 


Experimental Procedure 

An apparatus has been fabricated and tested 
which has been used to measure the thermal conduc- 
tance of pressed sample pairs at temperatures from 
1.6 K to 4.2 K under applied forces of up to 700 N. 
The apparatus is pictured in Figs. 1 and 2. In 
operation, it is immersed in a Dewar filled with 
liquid helium 4. To obtain data below 4.2 K, the 
temperature of the liquid helium 4 is reduced by 
evaporative cooling. A pressure controller limits 
evacuation to achieve the desired temperature. 

The following relation describes the mechanism 
of heat flow across the boundary between pressed 
solid surfaces: 

Q = k AT (1) 

where Q is the thermal energy transferred across 
the boundary, k is the effective thermal conduc- 
tance of the contact, and AT is the temperature 
difference across the boundary. While Eq. (1) is 
valid at any particular temperature T, k is 
actually a function of T so that: 

Q = J k(T)dT (2) 

If a simple power law is assumed to represent k(T) 
and temperatures Tu and Tl as measured at each 
side of the boundary, it follows that: 

/•Tu 

Q = J ctT n dT = a(Tu n+1 - TZ n+l )/(n + 1) (3) 

As discussed below a and n are determined empiri- 
cally for each sample pair, using Eq. (3). 

The actual pressed contact sample pairs were 
fabricated from OFHC copper. The samples were 
cleaned and stored in a nitrogen environment to 
prevent contamination while not under test. Five 
pairs were prepared to evaluate the effect of dif- 
ferent surface finishes on thermal conductance. 
Surface finishes of 0.1, 0.2, 0.4, 0.8, and 
1.6 urn rms were selected. Each sample pair had 
the following dimensions: 10.2 mm diam and 10.2 ram 
length for the lower sample, and 12.7 mm diam and 
8.89 ram length for the upper sample. The large 
diameter of the upper sample is to assure that any 
slight lateral movement would not prevent complete 
surface contact with the lower sample. Calibrated 
germanium resistance thermometers were installed 
in the upper and lower samples. A heater consist- 
ing of manganin wire wound on an aluminum form was 
placed above the upper sample (see Fig. 2). An 
analysis was performed to determine the losses of 
heater power due to radiation and gas conduction. 
This analysis showed such losses to be negligible, 
assuring that the measured heater power was the 
power actually applied. Each sample pair was 
tested from 1.6 to 4.2 K over the range from 0 to 


670 N. For each temperature and force value, 
heater powers from 0.100 mW to 10.0 mW were 
applied, the upper and lower sample temperatures 
measured. The force was then automatically incre- 
mented to the next value. In this manner, data 
sets were obtained for both ascending and descend- 
ing values of the applied force, to assist in 
determining the magnitude of a hysteresis effect 
(if any) . It should be noted that the apparatus 
was in all cases cooled down from room temperature 
with essentially zero applied force. Loading of 
the sample was performed during data acquisition 
only, and the load was relaxed when changing bath 
temperatures . 


Results 

Tables 1-8 present results obtained for the 
0.4 ym surface finish Cu-Cu sample pair. The first 
column in the tables denotes the bath temperature 
at which the data were taken. (Tables 1-7 list 
results for both ascending and descending force 
values, while Table 8 gives the ascending value 
only, since 670 N is the upper force limit.) The 
next two columns give the derived values of a 
and n. These values were obtained from a computer 
program which fit the experimental values of 
heater power Q, upper sample temperature Tu, and 
lower sample temperature TZ to Eq. (3) . The 
program also computes the magnitude of the offset 
in heater power Qo such that with no heater power 
applied a line fit of Q versus AT passes 
through the origin. The values of Qo are shown 
in column 3. In addition, the program performs a 
statistical analysis of the data in terms of the 
known uncertainties in the experimental measure- 
ments to calculate an uncertainty in the computed 
quantities a, n, and Qo. The next three columns 
give the standard deviation in the computed quan- 
tities a, n, and Qo as a measure of this error. 
The last two columns represent the minimum and 
maximum temperatures over which the computed values 
are accurate, representing the temperature range 
over which the original data were taken. 

In addition to the values of a, n, and Qo, 
Tables 1-8 provide the respective error associated 
with these values. By specifying the uncertainty 
in the experimental data due to measurement accu- 
racy and round-off errors, and perturbing these 
data assuming a Gaussian distribution, a standard 
deviation of the values 1 input to the computer 
program is obtained. Employing a random number 
generator, 99 computations of a, n, and Qo were 
performed by the program within the standard devi- 
ation in the output values of a, n, and Qo. In 
Fig. 3 the effective thermal conductance is plotted 
as a function of temperature for a range of applied 
forces. The curves shown in Fig. 3 were generated 
by averaging the values of a and n for both 
ascending and descending values of a particular 
force and plotting aT n as a function of 
temperature. 


Discussion 

The values of the exponent n given in 
Tables 1-8 correspond well with earlier work. 
Berman 1 observed a nearly T 2 temperature depen- 
dence of thermal conductance at liquid helium 
temperatures. The present range of values 
n = 1.9 to 2.2 certainly supports this finding. 
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Comparing the results in Fig. 3 with those of 
Berman 1 shows that in particular, at an applied 
force of 670 N (150 lb), the value of 8.0 x 10“ 3 W/K 
is within a factor of 2 of Berman's value 
1.46 x 10“ 2 W/K. Since Berman simply specifies a 
machine finish, it is felt that satisfactory agree- 
ment exists. Also shown is the effect of increas- 
ing applied force on thermal conductance. It is 
evident that the thermal conductance very definitely 
increases with increasing force, again supporting 
the earlier work of Berman. 1 


Conclusions 

It has been shown that for the 0.4 pm OFHC 
copper pressed contact pair, the thermal conductance 
varies roughly as the second power of the tempera- 
ture, and increases with increasing applied force. 

Future work will focus on copper sample pairs 
of differing finishes as well as stainless steel, 
aluminum, and brass and silica glass samples. Also, 
the effect of such coatings as gold and indium on 
thermal contact conductance will be investigated. 
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Table 1 Copper sample pair 0.4 um surface finish 
22 N applied force 


Table 2 Copper sample pair 0.4 ym surface finish 
45 N applied force 
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Table 3 Copper sample pair 0.4 pm surface finish 
112 N applied force 
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Table 6 Copper sample pair 0.4 pm surface finish 
448 N applied force 
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Table 4 Copper sample pair 0.4 pm surface finish 
224 N applied force 
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1.413 E-1 
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7.611 E-3 
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0.1564 

0.2167 
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1 084 E-1 
1 438 El 

6174 E-4 
1 090 E-3 

1773 E 3 
3 672 E-3 

5.936 E-4 
5.638 E-4 
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2 581 
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5245 

10 
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DESCENDING 

0.1550 

0.1966 

2.079 
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1477 E 2 
6473 E 2 

5.747 E-4 
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1653 E-3 
1846 E-3 

6 851 E-4 
6466 E-4 
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5685 

6200 

1.6 

ASCENDING 

OESCENOING 

0.1591 

0.1972 

2.122 

1128 

1.9940 E-2 
14153 E 2 

6 420 E-4 
6592 E4 

1509 E 3 
1576 E 3 

7 372 E-4 
7.043 E-4 

1 687 
1689 
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almW/K 1 ** 1 ! 

* 
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«°o 
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(KELVIN) 

3.8 
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0 3892 
0.3913 
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1603 E-2 
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1.320 E-2 
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4 858 
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3.6 

ASCENOING 

DESCENDING 

0 3758 
0.3939 
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1018 

8.104 E-2 
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9.743 E-3 
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1804 E-2 
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1198 E-3 
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14 
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12 
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0 3808 

1066 
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1447 E-1 
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1 937 E-3 
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11 

4532 

10 

ASCENOING 
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0.3529 
0 3596 

1079 

1081 

1 388 E-1 
1.402 E-1 
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4.821 E-3 

9.904 E-3 
1006 E-2 

1641 E-3 
1662 E-3 

2.997 

4 462 

14 

ASCENOING 

DESCENDING 

0 2798 
0.3029 

1063 

2.104 

1094 E -1 
2 .165 E l 

1604 E-3 
1 506 E-3 

4 400 E-3 
1910 E-3 

5.385 E-4 
4 992 E4 

2575 

2.406 

4 844 
4 683 

10 

ASCENOING 

OESCENOING 

0.7749 

1123 

•1.112 E-1 

1.064 E-3 

1096 E-3 

5 807 £4 

1 994 

4.700 

1.6 

ASCENOING 

OESCENOING 

0.2778 

0.2962 

1153 

1140 

8.2S19 E-2 
8808 E-2 

9.700 E-4 

1767 E-3 

6184 E-4 

1 689 

4528 


Table 5 


Copper sample pair 0 
336 N applied 


4 pm surface finish 
force 


Table 8 


Copper sample pair 0 
672 N applied 


4 pm surface finish 
force 
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4966 

3.6 
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7.727 E-2 
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3.4 
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DESCENDING 
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0 3388 

1.016 

2026 

1 328 E-1 
1 524 E l 

5.362 E-3 
6.265 E-3 

1.185 E-2 
1-306 E-2 

1.796 E-3 
1.952 E-3 
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3380 

4.821 

4.745 

3.2 

ASCENDING 

OESCENOING 

3077 

3287 

2.038 

2045 
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30 
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3.711 E-3 
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0.1860 

0.2566 

2078 
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1.292 E-3 
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2.3 
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•9460 £ 2 
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1.6 
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8.623 E 3 
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- 
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_ 

_ 
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2.151 
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2.794 E-3 
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Fig. 1 Overall view of thermal contact apparatus. 
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Fig. 2 Detail of cold plate. 
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Fig. 3 Results for 0.4 pm copper sample pair. 
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